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ABSTRACT 

We have developed a set of dynamically evolving Fokker-Planck models for the collapsed-core glob- 
ular star cluster M15, which directly address the issue of whether a central black hole is required to 
fit Hubble Space Telescope (HST) observations of the stellar spatial distribution and kinematics. As 
in our previous work reported by Dull et al., we find that a central black hole is not needed. Using 
local mass-function data from HST studies, we have also inferred the global initial stellar mass func- 
tion. As a consequence of extreme mass segregation, the local mass functions differs from the global 
mass function at every location. In addition to reproducing the observed mass functions, the models 
also provide good fits to the star-count and velocity-dispersion profiles, and to the millisecond pulsar 
accelerations. We address concerns about the large neutron star populations adopted in our previous 
Fokker-Planck models for M15. We find that good model fits can be obtained with as few as 1600 
neutron stars; this corresponds to a retention fraction of 5% of the initial population for our best fit 
initial mass function. The models contain a substantial population of massive white dwarfs, that range 
in mass up to 1.2 Mq. The combined contribution by the massive white dwarfs and neutron stars 
provides the gravitational potential needed to reproduce HST measurements of the central velocity 
dispersion profile. 

Subject headings: globular clusters: general - globular clusters: individual (M15) - stars: stellar dy- 
namics, mass function 
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1. INTRODUCTION 
1.1. Motivation 



As recently reviewed by iVesperiiil etall (IMoh . 
the possibility that globular clusters may contain 
intermediate-mass black holes (IMBHs) has received con- 
siderable attention, particularly in the past decade. An 
extrapolation of the relation between central black hole 
mass and velocity dispersion (Mbh — cr) for galaxies im- 
plies that globular clusters may harbor black holes with 
masses in the range of 10^ — 10^ M©. There have been 
a number of suggestions that the prototypical collapsed- 
core globular cluster M15 may contain a central black 
hole with mass of near 10'^ Mq. A central black hole has 
been invoked in Ml 5 to explain both the central sur face- 
brightness cusp (|Newell. Da Costa, fc Norrig [19761) and 
the central rise of the vel ocity dispersion profile observed 
both from the grou nd (iPeterson. Seitzer. fc CudworthI 
[19891) and by HST (IGerssen et al. I I20021I2003D . How- 



ever. iGrab horn et al. I (|1992| ). iDuU et al.l ()1997l ). and 
iBaumgardt (2003,) have presented dynamically evolv- 
ing models for M15 that provide good fits to the 
surface-brightness and velocity-dispersion profiles with- 
out the inclusion of a central black hole. Instead, 
in these models, centrally concentrated populations of 
massive white dwarfs and neutron stars dominate the 
gravitational potential in the central region. After 
some initi a l unce rtainty caused by plotting e r rors i n 
iDuU et al.l ([l997[ ). and rectified in iDuU et al.l (f20M ). 

[email: bmurphy@butler.edu| 



iGerssen et all (1200 31 found that the iDuU et al.l ()1997f ) 
mo dels do, in fact, fit the HST STIS velocities reported 
bv Ivan der Marel et al.l (|2002[ ). without the addition of 
a black hole. Using proper motions of stars iMcNamaral 
(j2003) also concluded there is little direct e vidence for an 
IMBH in M15. Using equilibrium models 'Chakr abartvl 
(i2006i) reac h ed a s imilar conclusion. iHo et al. (.20031) and 
iBash et all (|2008[ ) failed to detect X-ray or radio emis- 
sion coincident with the cluster center. This is not nec- 
essarily inconsistent with an IMBH, since the accretion 
rate and/or the accretion efficiency are likely to be low. 
Nonetheless, M15 does not appear to be a strong candi- 
date to harbor an IMBH. 

The population of ~ lO'' neutron stars included in the 
[Dull et all ()1997[ ) models follows from the adoption of a 
reasonable initial mass function (IMF), an inferred to- 
tal cluster mass of 4.9 x 10^ Mq, and the assumption 
that all of the neutron stars are ret ained after forma- 
tion. As we extensively discussed in iDull et al.l (|1997| ). 
this last assumption is the least certain, given the like- 
lihood that a substantial fraction of the neutron stars 
formed in cluster s will be immedia tely ejected by their 
large birth kicks. iPfahl et al.l (|2002[ ) have examined this 
issue, noting that while consideration of birth kicks sug- 
gests that the retention fraction should be in the range 
1 — 8%, X-ray evidence suggests that it may approach 
20% in some clusters, producing populations of up to 
1000 neutron stars in rich clusters. iDrukierl ()1996[ ) found 
similar estimates for retention fr a ctions of neutron stars 
born in binaries. IGerssen et al. I ()2003() argued that the 
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neutron star population adopted in the lDuU et alj (|1997[ ) 
models is implausibly large and thus that it is reasonable 
to effectively replace part of the inferred neutron star 
population with a central IMBH of mass ~ 2000 Mq. It 
is important to note, however, that our Fokker-Planck 
models only require that there be a population of nonlu- 
minous remnants with masses in excess of about 1.2 Mq; 
it is not necessary that all of these objects be neutron 
stars. 

The IMF plays a critical role in the dynamical evo- 
lution of a cluster. This is particularly the case for a 
core-collapsed cluster, in which extreme mass segrega- 
tion causes the massive white dwarfs and neutron stars 
to dominate the central region. The distribution of these 
massive remnants gives us a constraint on the upper end 
of the IMF. Dynamical models for clusters provide a key 
means of probing the distributions of these extremely 
faint stars, which presently elude optical detection. 

There is a significant increase in the quantity and qual- 
ity of th e data ava i lable for detailed modeling of M15 
since the lDull et al.l ()1997D models were developed. Thus, 
we have developed new Fokker-Planck models for M15, 
to fit these data, that address the issue of the necessity 
of a central black hole. 

1.2. Background 

HST imaging studies have provided a wealth of in- 
formation on globular cluster structure, stellar content, 
and kinematics. This is particularly true for the central 
regions of high density clusters, which were largely 
inaccessible from the ground before the HST-eia, owing 
to severe crowding of stellar images. This has allowed an 
examination of cluster m ass functions to unpre ce dente d 
depth "Sosm (199JbO; IParesce fc De Marchil ([200?} ): 
[Pasquali et al. (20Qj). In all three of these studies, 
multi-mass King models were fit to the mass-function 
and kinematic data to constrain cluster structure 
parameters. In less crowded outer regions of clusters, 
HST imaging has been used to detect stellar masses 
near the hydro gen burning limit of the ma i n sequence 
in N G C 6397 dPiotto. Cool, fc Kind [l997|: IKing et al ' 



a cluster (iMurphy fc Coh^ll988t iGrabhorn et 311 119921: 
iDuU et al 1 119971: jPasauah et al.l I2004D . More massive 



19981 [Richer et al.l l2008| ) and in M4 I Richer et al 



20021) . In the central regions of clusters, these improved 



observations allow us to investigate the details of 
the core collapse process. A predicted observational 
signature of core collapse is the presence of a central 
power-law surface densit y cusp surrounding a small, 
possi bly unresolved core (|CM [T98(H IMurnhv fc OAnl 
|1988[) . Ground-based observations of surface-brightness 
profiles have detected central cusps that are inter- 
preted as being due to core collapse in at least 21 
globu l ar clusters (see e.g. Jrag er. Di orgovski. fc Kingl 
Il99l ILugger. Cohn. fc Grindla^ [1993) ■ HST imag- 
ing studies of star count profiles have confirmed the 
presence of these cusps i n a num b er of clusters, includ- 
ing M1 5 (lYannv et all Il994al: iGuhathakurta et al. 
1991 ISosin fcKind I1997D. M30 (lYannv et al 



1994b; Sosin i997a)7 and NGC 6397 fSosinl 
Noyola & Gcbhardt (2006) have analyzed HST profiles 
for 38 clusters, finding that about half of the sample has 
a central power-law cusp, with a range of slopes. 

An important predicted outcome of the dynamical evo- 
lution leading to core collapse is that substantial mass 
segregation will occur within the half-mass radius of 



More massive 
stars will lose energy to lower mass stars and sink closer 
to the cluster center, as the stellar energy distribution 
tends towards equipartition. 

Indications that the stellar populations in the central 
cusp regions have been modified by close stellar encoun- 
ters provide other evidence of the occurrence of core col- 
lapse in clusters. Th ese include color gradients in the 
sense of inward bluing (jDjorgovski & Piotto''1993') , a cen- 
tral deficit of bright red giants (Stetson 1994), and cen- 
trally concent rated populations of bot h blue stragglers 
(iBailvn 1993) and millisecond pulsars ([AndersonI 119921 : 
iPhinnevllTggai . 

M15 (NGC 7078) has long been regarded as the pro- 
totypical core-collapsed cluster. In addition to its cen- 
tral surface-density cusp, it has other evidence of a 
high rate of stellar interactions in i ts central region in- 
cluding eight millisecond pulsars (lAndersonI I1992D. a 
blue straggler population ([de Marchi fc Parsecd 1199-^ , 
three low-mass X-ray binaries , AC 211, M15 X-2, and 
X-3 (Hertz & GrindlaiLl IT983I : IWhite fc AngelhH mM . 
iHein ke et al. 2009), and at least o ne cataclysm i c var i- 
able located in the cluster core IShara et al.| j200l . 
iDieball et all ([2001 and iHaurberg et al.l ([2O1O0 have 
found evidence for centrally concentrated populations of 
close binary candidates in Ml 5, including cataclysmic 

v ariables and helium whi te d warfs. 

iGrabhorn et al. I ([1992) and lDull et all ([1997) have fit- 
ted Fokker-Planck models to the surface-brightness and 
velocity-dispersion profiles of M15 and found that it was 
best fit with a global mass function having a power-law 
index oi x ^ 0.9, where 1.35 is the Salpeter index. These 
models also containned a sizable population of nonlumi- 
nous o bjects with masses in excess of 1 Mq. InLDuU^elalJ 
([1997f ) , the most massive of these objects were taken to 
b e neutron stars with mas ses of 1AM(^. 

iSosin fc Kiii^ ([1997D and lPiotto et al.l ([1997D found ev- 
idence for strong mass segregation in M15, using the 
HST Faint Object Camera (FOC) and Wide Field Plane- 
tary Camera 2 (WFPC2) imaging, respectively. At radii 
of 5", 20", and 4.5' they found that the mass func- 
tion had a power-law slope of —2.2, 0.2, 2.1, respec- 
tively, for stellar masses from 0.78 to 0.60Mq (where 
the Salpeter slope is 1.35). Significant mass segrega- 
tion has also been f ound in a number of other clus- 
ters, including M71 jRicher_fc Fahlm an 1989,) . 47 Tuc 
(Anderson fc King 1966; Howell et aLl lBol. NGC 6397 
Sosin. 1997bl and M30 (Sosin .1997aib[ ). ISosin fc Kingl 
19971) fitted a King-Michie model to the observed data 



set for M15 and found they could explain the observed 
mass segregation, but that they could not accurately re- 
produce the mass function in the centra l cusp nor the 
gravit ational potential of the cluster. iPasauali et al.l 
([20041) have supplemented the So sin fc Kingl ( 19971 ) mass- 
function dataset with HST NICMOS observations at 7', 
and fitted this combined dataset together with several 
kinematic datasets with a King-Michie mod el. The re- 
sulting fit is more succesful than that of Sosin fc Kingl 
(_1997i), most likely due to the inclusion of a much larger 
population of massive remnants. King-Michie models, 
which are non-evolving, may not accurately represent the 
current dynamical state of a post-collapse cluster such as 
M15. These models are likely to substantially underes- 
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timate the degree of mass segregation that occurs in the 
core and may thus provide a poor constraint on the rem- 
nant population. 

M15 provides us with an excellent opportunity to ex- 
amine the IMF of globular clusters. In all likelihood 
it has undergone little tidal str ipping due to i t s po- 
sition and orbit iii the Galaxy (iDauphole et al.l 119971 : 
iVesperini fc Heggid[l997l: iDinescu et al.lllQQQf ). The fact 
that M15 is a post-core collapse cluster also allows us to 
probe the upper main sequence of the IMF. Because the 
slope of the surface-brightness cusp is dependent upon 
the mass of the most massive stellar remnants we can 
estimate the number of n onluminous remnants present 
in M15 ()Murphv fc Cohn| [l988l. M15's low metallic- 
ity (Z=0.0001) implies that neutron stars were formed 
from progenitor stars of 6AMq and above (jHurlev et al.l 
l2000f l. From this we can estimate the slope of the IMF 
above this mass. The availability of detailed velocity pro- 
files and millisecond pulsar accelerations for M15 also 
aids in the estimate of the amount of mass in massive 
nonluminous re mnants ([Phinnevl 119931 : lAndersonI 119921 : 
IDull et al.llT99l . 

1.3. Overview 

In this paper we present dynamically evolving Fokker- 
Planck models of M15 that are more detailed than those 
of previous investigations. We fit these models to HST 
observations of stella r ma ss functions at seve ral radii by 
iSosin fc Kiilgl (jl997| ) and iPiotto et al.l (|1997D . In addi- 
tion we fit ou r mo d els to the v elocity dispersio n data 
of Gcbhardt etHI (Il99l [1997. 2000), DrukiereLaU 
([1998) , IGerssen et al. I (|2002[ ). and lMcNamaral (|2003[ ). In 
§2 we describe our models. The results and a newly de- 
termined best-fitting Fokker-Planck model for M15 are 
presented in §3. Section 4 summarizes our findings and 
discusses their implications. 

2. DYNAMICAL MODELS FOR M15 

2.1. Approach 

As in our previous m odeling of M15 (jGrabhorn et al. 1 
[1991 IDull etaJIfTOfll . our basic approach is to adopt a 
global stellar mass function, generate evolving Fokker- 
Planck models for this mass function, compare these 
models with the observed surface-density and velocity 
dispersion profiles at a number of epochs, and adjust the 
adopted mass function to achieve agreement with the 
observations. Our previous modeling of M15 adopted 
a power-law fo rm for most of the range of the stella r 
mass function (jGrabhorn et al. Ill992t iDull et al.lll997[ ). 
The direct observation of local stellar mass functions by 
HST studies provides an important new constraint in 
the fitting process, that allows the use of a more com- 
plex form for the mass function. Since HST observa- 
tions are now able to probe the stellar mass function to 
near the hy drogen burning limit of the main sequence 
(Ide Marchi & Parsccc 1994 : iCool. Piotto. fc Khlail996l: 
iKing et al.lll998t iRicher et al.ll2008| ). modeling these ob- 
servations requires good resolution of the mass function 
using 20 or more mass groups; in our previous model- 
ing we ha.d use d seven groups (jGrabhorn et al. 1 119921 : 
IDull et aI]|1997D . Our goal in the present study is to 
constrain the form of the global mass function by fitting 
observed local mass functions at several radii. 



Our globular cluster models are generated using the 
multi-mass direct Fokker-Planck method, which is based 
on a smooth, conti nuous description of the spatial struc- 
ture of a cluster (jMurphv fc CohnI 119881 ). In this ap- 
proach, each stellar mass group is represented by a sta- 
tistical distribution function in orbital energy. The dy- 
namical evolution of these distribution functions due to 
the effects of star-star gravitional scattering are followed 
in a diffusive approximation, in which the cumulative ef- 
fects of more frequent distant encounters dominate those 
of less frequent close encounters. The gravitational po- 
tential due to the stellar distribution is determined self- 
consistently, yielding the time evolution of the structure 
of the model globular cluster. Using this method we 
were able to follow the evolution of a set of initial mod- 
els through the first core collapse and the subsequent 
post-collapse evolutionary phase using heating via dy- 
namically formed bin aries (iMurphv et al.lll990l). As re- 
viewed bv lGoodma^ (11993D andjMurphv (199f), energy 
input into the cluster core from interactions involving 
close binaries is expected to halt core collapse and drive a 
post-collapse oscillatory phase. Binaries may be present 
in the cluster from primordial formation and/or may be 
formed by two-body tidal capture and three-body dy- 
namical processes. M15 likely has few binaries in its core 
due to its being in deep collapse. Deep collapse of a 
cluster is an indication that it has "burned" most o f the 
primordial binaries near its center (IGao et al.lH99l . 

We begin our models with a mass-density distribution 
represented by a Plummer model (n = 5 polytrope). 
All models start with the mass groups having identi- 
cal vel ocity distributions, a,s expe cted after violent relax- 
ation (jBinnev fc Tremaind IT9871 §4.7.2). This velocity 
equipartition is quickly replaced by energy equipartition, 
causing the most massive stars to sink to the cluster cen- 
ter as the cluster evolves dyn amically (jlnagaki fc Saslawl 
119851 IMurphv fc Cohnll 19881 ) . Large-scale central density 
oscillations appear soon after the initial collapse. These 
oscillations have been observed in conducting gas-sphere 
models, Fokker-Planck models, and iV-body simulations 
jSugimoto fc Bettwieser"1983l: IGohn. Hut, fc Wis(j[l989l : 
[Makino. ,1996) . As discussed in our previous work, the 
core radius lies close to its maximally expanded size over 
most of an oscillation cycle. For a cluster at a distance 
of ~ 10 kpc, such as M15, this corresponds to an angu- 
lar scale of ^ 1 — 2", which is similar to observational 
upper li mits on the core radi us of M15 from HST imag- 
ing (e.g. ISosin fc Kingl 119971 ). We typically evolve our 
models through core collapse and into the post-collapse 
phase, to provide a wide choice of evolutionary phases 
with which to fit the observations. 

Recently Fregeau, (2008) has suggested that clusters 
like M15 may still be in a binary burning state rather one 
of deep collapse; we note that this results in a similar core 
radius to that achieved in a post-collapse core bounce. 

2.2. The Stellar Mass Function 

Our models of M15 approximate the present-day stel- 
lar mass function with a set of 20 discrete evolved mass 
groups. This mass resolution has allowed us to inves- 
tigate the effects of mass segregation on the local mass 
function. Rather than explicitly include continuous stel- 
lar evolution mass loss, which is most important in the 
early evolutionary phases, we use a static present-day 
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TABLE 1 

Characteristics of the Stellar Mass Functions 



Bin 


Stellar 


Mass 


Type 


Progenitor 




Mass (M0) 


Fraction 




Mass {A/0) 


1 
1 


1 y1 /I n 


0. 


.0051 


neutron star 


P. A'x 00 n 


9 


1 91 fl 


0. 


.0684 


white dwarf 


A nn +r» R A'\ 


o 
O 


u.you 


0. 


.1339 


white dwarf 


nn -f^-. A nn 

Z.UU TO 


A 


U. I UU 


Q 


.1840 


white dwarf 


1 nn nn 
i.UU TO z.uu 


U 


U.u /U 


0. 


.0638 


white dwarf 


n CO \^r^ 1 nn 
u.oz to l.UU 





U.oiD 


0, 


.0036 


red giant 




/ 


U. / yD 


0, 


.0209 


main sequence 




Q 
O 


n VQn 
U. / oU 


0, 


.0674 


main sequence 




g 


0.619 


0, 


.0963 


main Gf^mi/^nfp 




10 


0.516 


0, 


.0620 


main sequence 




11 


0.442 


0, 


.0420 


main sequence 




12 


0.378 


0, 


.0326 


main sequence 




13 


0.323 


0, 


.0320 


main sequence 




14 


0.276 


0, 


.0315 


main sequence 




15 


0.236 


0, 


.0310 


main sequence 




16 


0.202 


0, 


.0306 


main sequence 




17 


0.173 


0, 


.0301 


main sequence 




18 


0.148 


0, 


.0281 


main sequence 




19 


0.126 


0, 


.0211 


main sequence 




20 


0.108 


0, 


.0154 


main sequence 





mass function for all epochs of our models, derived from 
the IMF by transforming post-main-sequence stars to 
appropriate remnant masses. This present-day evolved 
mass function takes into account stellar evolution that 
has modified the upper IMF over the lifetime of the clus- 
ter. For this study we were mainly interested in the final 
remnant mass of stars so that we could determine the 
number and mass of the objects in each of our remnant 
mass bins. To determine the initial main sequence to fi- 
nal remnant masses for these bins we used t he ste llar evo- 
lution models described by iHurlev et al.l (|200Clf ) . From 
their fitting formula we were able to determine the mass 
of the remnants produced by upper main sequence stars 
of the IMF. Five of the twenty mass groups represent 
degenerate stars, two represent post-main sequence stars 
(horizontal branch and red giant), and 13 mass groups 
represent the main sequence stars over the mass range 
0.8 to O.lMrn Usi ng a metalhcity of Z^O.OOOl in the 
IHurlev et al.l ()200C1[ ) models results in significantly more 
massive white dwarfs than would solar metallicity. Also 
neutron stars are produced from progenitor stars as low 
as 6.43Mq from this low metallicity. 

The IMF of the cluster is given by a series of power- 
laws of the form, 

(i+^i)dm, 22M0 > TO > mi; 
m^^^+^^^dm, TOi > TO > TO2; 

-(l+a:3)(^^^ m2 > TO > TO3 ; 



dN 



TO 



(1) 



^ ■) ^ ' ' o; 

-(1+3:4)^^^ m4 > m > O.IMq 



Here dN is the number of stars contained in the mass 
interval {m,ra + dm); xi, X2, X'^, and Xa are the loga- 
rithmic slopes of the mass function over the indicated 
mass intervals Q; and toi, TO2, TO3, and TO4 are the stel- 
lar masses where the initial mass function is allowed to 
change power-law slope. For progenitor masses above the 
turnoff mass, we assumed a single power-law slope that 
was varied from x = 0.9 to 1.3, where the Salpeter mass 
function is a; = 1.35. To estimate the number of stars in 

^ The parameter x is defined as a: = —dlnNjdlnm. 



the five degenerate mass groups, we integrated the IMF 
over the appropriate progenitor mass ranges as indicated 
in Table 1. As the most massive group of degenerate 
stars, neutron stars of mass 1.44AfQ will determine the 
slope of the surface-density profile and the local mass 
function in the inner parsec of the cluster if they domi- 
nate the mass density in the cusp. We assume that neu- 
tron stars originate from stars with initial masses from 
6.43 t o 22 solar masses as prescribed by IHurlev et al.l 
(l2.00Cf) . The properties of each stellar group are listed in 
Table 1. 

Changing the upper mass limit for neutron star pro- 
genitors has only a small effect on the total number of 
neutron stars produced. Becaus e neutron stars receive 
a sizeable velocity kick a t birth ()Lvne fc Lorimeil 119941 : 
iHansen fc Phinnev 1 119971 ) . it is likely that a significant 
fraction has been ejected from the cluster. We there- 
fore allowed the number of neutron stars to be lower 
than would otherwise be presumed from the IMF. For 
a globular cluster similar to M15, the central escape ve- 
locity is typically on the o rder of 50 kms~^ . For t he pul- 
sar velocity distribution of lLvne &: Lorimeil (|1994l ). fewer 
than 10% of the o riginal neutron stars may be present 
(jPfahl et all l2002| ). This quantity might be higher if 
one assumes tha t some of the neutron stars orig i nated 
in binary stars (iDrukieii 119961 : iDavies k, Hansen! 119971 : 
iPfahl et al.l l2002f) . An important goal of the present 
study has been to determine the minimum retention frac- 
tion for which we are able to fit the M15 observations. 
For other degenerate groups we assumed that all stars 
are retained in the cluster. The stellar mass of each 
white dwarf bin is given in Table 1, along with the as- 
sumed progenitor masses. The mass of t he white dwarf 
bins were selected using a prescription of iGirardi et al.l 
(2000); Hurley ct al. (2001). Using this method, we de- 
termined the number of remnant stars present by inte- 
grating the IMF over the various range of masses above 
the current turnoff mass. 

Stellar mass black holes are not included in these mod- 
els for several reasons. As with neutron stars bl ack holes 
may receive a birth kick (Dhawan "etall[2007l ). If this 
does occur then very few black holes would be retained 
by the cluster. Even if there were no birth kick most 
of the black holes produce by a typical mass function 
would likely be ejected. Given the mass difference com- 
pared to the rest of the cluster stars these stellar mass 
black holes would undergo rapid dynamical evolution, 
quickly segregate to t he cluster core, and decouple from 
the rest of the cluster (|Spitzeiill969HWatters eraI]|2000D . 
This central cluster of black holes would then form bi- 
naries that would then cause the m to be eje cted from 
the system a relatively short time (jKulkarni et al... 19931 : 
ISigurdsson fc Hernquistl[l993l ). 

2.3. Model Fitting 

Our models can be described by three different param- 
eters, the present-day total mass of the cluster, the initial 
half-mass radius of the cluster, and the stellar mass func- 
tion parameter se t. We first considered an IMF similar 
to that inferred bv lDull et all (|1997f ). The Fokker-Planck 
model was first run with this mass function, then the re- 
sults were examined to determine if a fit could be made, 
at some epoch, to all three local mass functions found by 
ISosin fc King. (,1997. ) and .Piotto et al. (,1997. ) at 5", 20", 
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1 0.8 0.6 0.4 0.2 0.1 



Mass Mq 

Fig. 1. — The fit to the stellar mass function of M15. From top to 
bottom the radius of sampling for the mass functio ns are 5", 20", 
and 4 .5 '. Th e observational data arc adapted from fSosin fc Kind 
119971 ): [Fiott o ct al. (1997). The model fit at 4.5'is an average of 
three radii, 3.90', 4.28', and 5.46'corresponding to HST WFPC2 
chips W2, W3, and W4, respectively. 

and 4.5'. From this first model we were able to gauge 
what adjustments were necessary to the IMF, the half- 
mass radius, and the cluster mass. We revised the model 
until it satisfactorily matched the mass-function obser- 
vations at the three radii at some epoch. The models 
were fit to four types of data: the mass-function data at 
the three radii, star-count profiles for three luminosity 
ranges, the velocity-dispersion profile, and millisecond 
pulsar acceleration measurements. A model was deemed 
acceptable if it satisfactorily matched all four of these 

data sets. 

Sos in fc Kind ()1997[) present radial star-count profiles 
for several magnitude ranges. For a comparison with 
their results we converted our stellar mass functions to lu- 
minosity functions using the mass to lumi nosity relation 
(MLR) of lBergbusch fc Vandenbeii (1199^ . We assumed 
a dist ance modulus of 15.4 identi cal to iSosin fc Kind 
(|1997f) and lDurreU fc Harris I (|1993[ ). These MLRs allow 
us to make a direct co mparison to the star count data 
of ISosin fc Km^ (fT997[) at radii of 5" and 20 ". The 
MLR is uncertain for low metallicities and low masses 
([Alexander et al.lll997D . so we concentrated on fitting the 
mass function for stellar masses above 0.4Mq. 

3. MODELS 

The best-fitting model has a mass of 4.4 x IQ^AIq and 
a half-mass radiu s of 3.4 parsecs. This m ass is similar 
to that found bv iMcNamara et al.l ()2004D . The model 
was found to fit best in the post-core collapse phase. 
While a fit of the model at the initial collapse gave rea- 
sonable results, we found that the post-collapse phase 
provided a better fit to both the local mass functions 
and the star-count data. We tried to fit a number of 
models at several points during a post-collapse phase, in- 
cluding both expansion and recollapse phases. We found 
that the maximally expanded model could not account 
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Fig. 2. — Star count p rofiles for three luminosity ranges compared 
to observational data of lSosin & Kinel llT997'). The lower solid curve 
is for V magnitudes between 18.25 and 19.75 (M = 0.80 Mq). The 
upper solid curve is for V magnitudes between 19.75 and 21.25 
(M = 0.73 Mq). The dashed curve is for V magnitudes between 
21.25 and 22.75 (M = 0.62 Mq). 



TABLE 2 
M15 Initial Mass Function 



Power-law Slope 


Mass Interval 


X 


(Mq) 


0.9 


m > 0.55 


-1.5 


0.55 > m > 0.40 


0.9 


0.40 > m> 0.15 


-1.0 


0.15 > m > 0.10 



for the measured pulsar accelerations. In addition to 
the finding the best epoch to fit the model we found the 
slope of the upper IMF is constrained by the number 
of neutrons and massive white dwarfs needed to fiatten 
the surface-brightness profile of the cluster and produce 
a large enough pulsar acceleration. The IMF slope could 
be steepened but then an unrealistic percentage of neu- 
tron stars would have to be retained. The break points in 
the mass function below the turn-off mass are primarly 
constrained by the slope of the observed mass function at 
4.5'. Break points significantly different than those listed 
in Table 2 led to a poor fit of the observed mass function 
at all radii. Our best-fitting model, which corresponds 
to an expanding phase, is shown in Figures 1 — 7. The 
mass-function slopes used over the four mass intervals, 
described in section 2.2, are given in Table 2. Above 
O.55M0 we assume a slope of x = 0.9 with a neutron 
star retention of 5%. The 1.4Mo objects in the best- 
fitting model make up 0.5% of the cluster mass. In the 
range of O.55M0 to Q.AMq the mass function flattens 
significantly, to x = —1.5. This trend reverses for the 
lower main- sequence, below O.4M0, where x increases 
once again to 0.9. Finally, below O.15M0, the mass func- 
tion flattens once again to a; = — 1. 



6 



Murphy et al. 



TABLE 3 

M15 Best-Fitting Model Parameters 



Parameter 



Value 



Total Mass 

Tc 

rh 

No. of Neutron Stars 

Projected Half-Mass Radius of Neutron Stars 
No. of Massive White Dwarfs (m = 1.21Mq) 
(m,) 

% Cluster Mass in Remnants 



4.4 X 10^5 Af0 
0.027 pc 
3.4 pc 
0.008 
1560 
0.17 pc 
24872 
0.43 Mq 
45.5% 



Figure 1 shows the local m ass function at radii of 5", 
20", and 4.5'. Note that the iSosin fc Kind ([1997) data 
only show stellar masses down to 0.6 and 0.5Mq for radii 
of 5", 20", respectively, due to observational limitations 
in the cluster cusp described earlier. Our model accu- 
rately reproduces the mass function at all three radii. 
Figure 2 shows star-count profile of the mode l com - 
pared to the FOG observations of ISosin fc Kingl (|1997f l 
for three magnitude intervals. The three magnitudes 
intervals correspond to subgiants and stars near the 
turnoff (18.25 < V < 19.75) and two main-sequence 
groups (19.75 <V < 21.25 and 21.25 <Y < 22.75). 
The model once again accurately reproduces both the 
slope and normalization of the profile. F igure 3 shows 
the v el ocity-dispersion measu remen ts of iDrukier et alJ 
(fT998h . IGerssen et ahl (pOOl . and IMcNamaral (|2008n 
compared to the model. The model passes through most 
of the error bars of the velocity p oints. We do n ot try 
to fit the outer flattening of the IDrukier et al.l ()1998f ) 
dataset, because this is likely to be an effect of being near 
the tid al radius. The King-Michie model of Sosin 
(|1997| ) matches the surface density well but provides only 
a fair fit to the mass function, and fails to reproduce 
the velocity dispersion profile. The King-Michie model 
of ([Pasauali et al.|[2004D matches the surface brightness 
profile and the luminosity functions, while also provid- 
ing a reasonable fit to the velocity dispersion profile. As 
discussed below, this appears to be a consequence of re- 
taining all of the neutron stars produced from the IMF. 
Our model has a much more realistic retention fraction 
of 5%. This reduced requirement for heavy remnants, 
to provide the central potential, is a consequence of the 
higher degree of central concentration that develops in 
our evolving Fokker-Planck models relative to a King- 
Michie model. Figure 5 shows the cumulative projected 
mass as a function of radius for the total mass and sev- 
eral of the mass groups. Using observed kinematic data 
for M15 Chakrabarty (2006) estimated a mass interior 
to 0.010 pc of roughly lOOOAf©. We find somewhat 
less mass interior to thi s radius in our model. How- 
ever iChakrabartvl ()2006D had an inferred total cluster 
m ass nearly 3 times that of our best- fitting model and 
of iMcNamara et all ()2004D . This discrepancy could be a 
result of mass segregation and our use of evolving models 
and the presence of heavy remnants. 

Eight milliseco nd pulsars have been detected in M15 
(|AndersonlfT99^ . Two of these, PSR 2127-l-llA and 
PSR 2127-l-llD, are located 1" from the cluster center 
and have measured accelerations. These accelerations 
allow us to determine a lower limit on the gravitational 
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Fig. 3. — Velocity dispersio n of best- f itting model of M15 
compa red t o observation s of Dr ukier et al.l II1998I ) , IGerssen et al. I 
ll200a) . and IMcNamaral 112003 '). Note that the central concentra- 
tion of massive white dwarfs and neutron stars produces a central 
one-dimensional velocity of 13.6 kms^^. 
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Fig. 4. — Mass surface density profiles for the best fitting M15 
model. Note that the region within 1' of the cluster center is dom- 
inated by white dwarfs, with massive {M > 0.95 Mq) white dwarfs 
dominant within the inner few ". Neutron stars are not dominant 
at any radius, due to the relatively modest retention fraction of 
5%. 



Non- 



potential (jPhinnevI [l99l IGrabhorn et al. 1 119921) . 
luminous remnants dominate at this radius, as seen in 
Figure 4. These remnants (all white dwarfs plus neu- 
tron stars) make up 45% of the total cluster mass glob- 
ally, and dominate the inner regions due to their greater 
mean mass than that of the luminous stars. We use the 
measured pulsar accelarations to determine the minimum 
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Fig. 5. — Cumulative projected mass for the best fitting 
model. The projected half-mass radii of the neutron stars 
total cumulative mass are 0.056' and 0.83', respectively 



10' 



M15 

and 



mass of nonluminous remnants required in our model. 
Any dynamical model must fit these accelerations in or- 
der to properly map the potential in the cusp. Figure 
6 shows the projected gravitational acceleration versus 
radius for the best-fitting model. A m aximally expande d 
model does not fit the accelerations ()Dull et al.lll997D . 
The pulsar accelerations combined with the fit to the 
velocity-dispersion profile gives us confidence that we 
were properly estimating the nonlumi nous remna n t pop - 
ulation within the half-mass radius. iDull et alJ (pOol 
found that for high-concentration King Michie models, 
the gravitational acceleration falls s hort in the inn e rmost 
regions of the cluster. The model oflSosin fc Kind (|1997[ ) 
(labeled "King" in Fig. 9 of DulLeJ^al. 2003) failed to fit 
the accelerations (as well as the central the velocity dis- 
persion profile) because they had very few massive rem- 
nants in their model, 21% at O.SSM© and less than 1% in 
a 1.55Mq bin. Because most of their remnants were in 
a lower mass bin, very few nonluminous r emnants segre- 
gated into the cusp region of the cluster. iPasquali et alJ 
|20G4) did not consider the pulsar accelerations. 

Because mass segregation modifies the observed mass 
function, one would be curious to see whether there is any 
region of the cluster where the local mass function is a 
close approximation to the global mass function. Fig- 
ure 7 shows the variation in the observed mass func- 
tion with radius. The observed mass function changes 
from being much flatter than the global mass function 
in the inner part of the cluster to being steeper than 
the global mass function in the outer regions of the clus- 
ter, as would be expected. None of the curves exactly 
rep licate the glob a l mas s function. However, as noted 
bv IPasauah et al.l (|20Q4D . the local mass function best 
matches the global mass function near the half-mass ra- 
dius. 

ft should be pointed out that at the outermost radii 
of our model cluster, i.e., radii greater than 10', the lo- 
cal mass function and global mass function are identical. 
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Fig. 6. — The line-of-sight acceleration as produced by the mean- 
field gravitational p otential of the b est-fitting model compared to 
the observations of lAndersori (1992). The data symbols indicate 
the minimum line-of-sight acceleration needed to account for the 
negative period derivatives of PSR 2127-l-llA (lower point) and 
PSR 2127-l-llD (upper point); the pulsars lie at nearly the same 
distance from the cluster center. 

This is due to the very long relaxation times at these 
large radii. We do not show the mass function here be- 
cause our cluster is assumed to be isolated whereas for a 
more realistic model some modification would likely have 
occurred at the outermost radii due to tidal stripping by 
the Galaxy. 

4. DISCUSSION 

Our inferred IMF is flatter than a Salpeter IMF over 
the entire initial mass range, with a slope of x = 0.9 for 
M > 0.55 Mq. We were able to fit the star count pro- 
files, local mass functions, kinematic data, and millisec- 
ond pulsar accelerations with a neutron retention frac- 
tion of 5% for neutron stars created for stars with ini- 
tial masses above 6.43M0. Without objects with masses 
of 1.2M0 and above, the star-count profiles would have 
been too steep over the radial range of to 10". The is- 
lope of stellar density profile of the most luminous stars 
depends on the ratio of mass es of the lum inous stars and 
the dominant mass species (fCoiml 119851) ■ This IS given 

by 



d IncTu, 

a = -— ^ = -1.89 

d mr 



mdo 



0.65 



(2) 



Where a is the logarithmic slope of the projected density 
profile, is the mass of the luminous stars (O.8IM0), 
and mdom is the mass of the dominant mass species 
(1.21Mq). Note that this equation implies that the slope 
of the cusp giants should be — 0.63, which is close to the 
obser ved slope (Sosin &: King|[T99 7: Noyola & Gebhard^ 
I2006f ). This slope applies to the region where massive 
remnants (M > I.2M0) dominate the surface density. 
This corresponds to a radius inside of 5" for our models, 
as seen in the surface-density profiles of our best-fitting 
model shown in Figure 4. Outside of this radius the 
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Fig. 7. — Radial variation of the stellar mass funetion as com- 
pared to the global mass function (dashed curve) . Note global and 
local mass functions are similar only near the half-mass radius. In 
this model = 3.4 parsec. 

surface-density slope slowly steepens and does not have 
a power-law form. 

The models presented here require a total population 
of about 1600 neutron stars, which is a dr amatic re- 
ducion relative to the 10^ neutron stars in the lDuU et alj 
(|1997t ) models. This appears to be much more real- 
istic. Nonetheless, it is possible that some of the 1.2 
and 1.4Mq objects may be represented by nonluminous, 
tightly bound white-dwarf binaries. In our models 32% 
of the total mass of the cluster is in white dwarf bins 
with masses 0.7 and 0.95 Mq. It could be that these 
objects have undergone exchange reactions with primor- 
dial binaries composed originally of lower-mass main- 
sequence stars, via 3-body or binary-binary interactions 
(jSigurdsson & Phinney 1993). This would not lower the 
number of single white dwarfs in the core significantly, 
due to their large numbers, but would add to the 1.4Mq 
bin. These binaries could not be present in large num- 
bers in the form of main-sequence pairs. Such main- 
sequence binaries would have an magnitude of roughly 
19 and would be apparent in the star count profiles of 
ISosin fc Kinal (l997). for the magnitude range 18.25 to 
19.75, in the form of a steeper profile. If main-sequence 
binaries more massive than the turnoff did represent a 
sizeable fraction of the stars in this magnitude range, 
they would cause the profile to be much steeper than is 
observed. 

One last explanation for the central slope of the sur- 
face density profile of Ml 5 is the possible presence of 
a massive central black hole. The slope of the star 
count profile tends to be near —0.7. This is very sim- 
ilar to the —0.7 5 exp ected from a massive black hole 
(jBahcah fc Wolj[l976h . In the muhi -mass case the slope 
of the projected star-count profiles in a Keplerian poten- 
tial are given by 

/ m, 1\ 
V4mrfom 2y 



(|Murphv. Cohn. fc Durisenlll991[ ). Note that if objects 
more massive than the turnoff stars dominate the stellar 
density, then the turnoff stars wih hav e a slope slightly 
less than -0.75 (jBahcah fc Woljll977[ ). Also note that 
the minimum slope possible for any group cannot be less 
than —0.5 if a m assive black hole is r esponsible for the 
cusp. Though the ISosin fc Kind (jl997[) data do show this 
trend, it should be pointed out that in order to achieve 
these slopes out to 10", the central black hole would 
have to have a mass in excess of 5OOO7\f0. Other evi- 
dence against a massive central black h ole comes from 
the ve locity -dispersion n ieasur ements of iGerssen et al. I 
(|2002f) and iMcNamaral (j2003D . which show a gradual 
rather than a steep rise in the velocity dispersion. The 
velocity dispersion in the Keplerian potential of a mas- 
sive black hole is expected to rise as ti oc r~°-^ as op- 
posed to w oc r~°-^. Thus, there is little need to in- 
voke the presence of an IMBH in M15 on the basis of 
ou r model fitting. This i s in a ccord with the conclusion 
by Ivan den Bosch et al.l (|2006D that the significant cen- 
tral dark mass concentration in M15 can either be in the 
form of an IMBH, individual objects, or a combination 
of these. 

One useful aspect of our models is that they can help 
guide obsevers as to what radii should be observed in 
order to best retrieve the global mass function. We find 
that observations near the half-mass radius give the best 
results for two reasons. First mass segregation has not 
altered the mass function as much as in the inner regions 
and second, tidal effects have likely had less effect at this 
radius than further out. Obtaining data at center of the 
cluster, 1 parsec from the center, and at the half-mass 
radius appear to be the best method for retrieving the 
global IMF. These radii are nearly identical to the data 
we have fit from Sosin & King (1997)- For clusters that 
are nearer to the Sun it is irn p ortant that observations 
be spaced accordingly. iSosinI (|1997b[ ) shows an exam- 
ple of what may happen when local mass functions are 
determined at two adjacent radial bins in the cusp of 
a cluster. For M30 (NGC 7099), Sosin examined local 
mass functions in bins centered at 3" and 9" from the 
cluster center, which corresponds to 0.1 and 0.3 parsecs, 
respectively. These radii are similar to our second and 
third curves in Figure 7. Figure 7 shows that the slope of 
the mass function above O.4M0 remains nearly constant 
at these small radii, even the though relative numbers of 
stars have changed. This should not be looked at as a 
sign of little mass segregation. This same effect would be 
seen even for a cluster in deep core collapse or bounce. 
This highlights the importan ce of radially s pacing obser- 
vations in globular clusters. iSosinI (jl997bf ) did obtain a 
HST FOG field at 20", but unfortunately the field was 
not usable due to the presence of a brigh t giant. 

Comparing our results with those of iPasquali et al.l 
(2004j), we obtain a similar slope of the IMF above 
0.55Mq. Below this mass, there are detailed differences 
in the mass function, but an overall trend of flattening or 
turnover. For both sets of models, there is a substantial 
amount of mass in remnants. For our models, the total 
mass in remnants is 46% of the cluster mass with 21% 
of the cluster mass in remnants above the turnoff mass. 
Given the evolving nature of our models, they are likely 
to provide a more realistic description of mass segrega- 
tion and the central gravitational potential of M15. 
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5. CONCLUSIONS 

We have fitted dynamically evolving Fokker-Planck 
models to a wide variety of data for M15, including HST 
observations of the stellar mass function at several radii, 
the star-count profiles, the velocity dispersion profile, 
and radio measurements of millisecond pulsar acceler- 
ations. We infer that M15 has a relatively flat IMF slope 
of 0.9 or flatter. With a modest retention fraction of 
5%, we infer that M15 presently contains 1560 neutron 
stars. While of order this amount of IAMq objects is re- 
quired to fit the data, some of these may be in the form 
of tightly-bound white dwarf pairs. The central cusp re- 
gion is dominated by 1.2Mq white dwarfs, which repre- 
sent about 7% of the cluster mass. For a low-metallicity 
cluster like M15, a substantia l number of white d warfs 
are formed in this mass range iHurlev et all (pOOOl ). Be- 
cause of mass segregation, the global mass function is not 



replicated at any point in the cluster. However, the lo- 
cal and global mass functions come closest to each other 
near the half-mass radius of the cluster. 

We conclude that there is no need to invoke the pres- 
ence of an IMBH in Ml 5. The presence of massive white 
dwarfs and neutron stars provides the necessary central 
gravitational potential to fit the surface density profile, 
the velocity dispersion profile, and the millisecond pulsar 
accelerations. 
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